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I Abstract 
> 

O We report the discovery of CoRoT 102980178 (a = 06'^50'"12!10, S = -02°41'21'.' 
J2000) an Algol-type eclipsing binary system with a pulsating component (oEA). 
It was identified using a publicly available 55 day long monochromatic lightcurve 
from the CoRoT initial run dataset (exoplanet field). Eleven consecutive l'."26 
I deep total primary and the equal number of 0™25 deep secondary eclipses (at 
phase 0.50) were observed. The following light elements for the primary eclipse 
were derived: HJDMini = 2454139.0680 -h 5'?0548 x £ The lightcurve modeling 
leads to a semidetached configuration with the photometric mass ratio q = 0.2 
and orbital inclination i = 85°. The out-of-eclipse lightcurve shows ellipsoidal 
variability and positive O'Connell effect as well as clear 0'."01 pulsations with 
5— i the dominating frequency of 2.75 c/d. The pulsations disappear during the pri- 
mary eclipses, which indicates the primary (more massive) component to be the 
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* The CoRoT space mission was developed and is operated by the French space agency 
CNES, with participation of ESA's RSSD and Science Programs, Austria, Belgium, 
Brazil, Germany and Spain. 
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pulsating star. Careful frequency analysis reveals the second independent pul- 
sation frequency of 0.21 c/d and numerous combinations of these frequencies 
with the binary orbital frequency and its harmonics. On the basis of the CoRoT 
lightcurve and ground based multicolor photometry, we favor classification of 
the pulsating component as a 7 Doradus type variable, however, classification 
as an SPB star cannot be excluded. 

Individual Objects: CoRoT 102980178 

Introduction 

Pulsating stars in eclipsing binary systems are objects of considerable astro- 
physical interest. It is well known that in double-line binaries masses, radii and 
luminosities of the components may be directly determined which allows one 
to constrain pulsation models. Possible influence of one component of a binary 
system on pulsations of the other component (through mass transfer or tidal 
interaction) is an interesting and still largely unexplored topic. A case of tidal ex- 
citation of 7 Doradus type pulsations was reported by Handler et al. (2002; see 
also references to some theoretical discussions therein). Finally, binary systems 
with oscillating component offer an intriguing possibility of direct observations 
of radius change during the pulsation cycle as periodically changing moments 
of 1st, 2nd, 3rd and 4th contact (in the case of radial pulsations), and in gen- 
eral case - possibility of pulsation mode identification by use of the secondary 
component as a spatial filter (Biro &i Nuspl 2005). Altogether, identification 
of pulsating stars of different types in eclipsing systems is important. 

The combination of pulsations and eclipses is often found among the evolved 
stars like symbiotic binaries (e.g. Chochol &i Pribulla 2000) and cataclysmic 
variables (e.g. Araujo-Betancor et al. 2005). Four Type I or II Cepheids in 
eclipsing binaries have recently been identified (see Antipin et al. 2007 and 
references therein), the confirmation of one more possible object of this type 
is still pending (Khruslov 2008). The discovery of a red giant with solar-like 
oscillations in a long-period eclipsing binary system was recently announced by 
Hekker et al. (2010). 

Semi-detached eclipsing binaries with a pulsating primary component which 
is close to the main sequence were dubbed "oscillating Algols" (oEA, Mkrtichian 
et al. 2002). About twenty such systems are known to date (Mkrtichian et al. 
2007) most of them containing 5 Scuti type components. The first eclipsing 
binary containing 7 Doradus type pulsating star (VZ CVn) was reported by 
Ibanoglu et al. (2007). Recently Damiani et al. (2010) and Maceroni et al. 
(2010) identified two 7 Doradus candidates in eclipsing binaries using CoRoT 
photometry (CoRoT 102931335 and CoRoT 102918586). 
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CoRoT (Convection, Rotation and planetary Transit; Fridlund et a! 2006) 
is a space experiment devoted to study of convection and rotation of stars 
and to detection of planetary transits. CoRoT data become publicly available 
one year after release to the Co-Is of the mission from the CoRoT archive: 
http: //idoc-corot . ias.u-psud.fr/ 

In this paper we announce the discovery of a new oEA system identified 
using publicly available data from the CoRoT initial run (exoplanet field). An 
almost uninterrupted 55 day long monochromatic CoRoT lightcurve together 
with ground-based multicolor photometry allow us to identify the pulsating 
component as a likely 7 Dor type variable, however classification as a Slowly 
Pulsating B (SPB) star can not be ruled out. The object is an interesting target 
for a spectroscopic follow-up. 

Observational data 
CoRoT photometry 

CoRoT 102980178 (USNO-Bl.O 0873-0161681, coordinates: a = 06^50™12!10, 
5 = -02°41'21'.'8, J2000, Monet et al. 2003) was identified by us as an Algol 
type eclipsing binary with oscillating component (oEA star) after visual inspec- 
tion of the CoRoT lightcurve. The star was also independently identified as an 
eclipsing binary by Debosscher et al. (2009) and Carpano et al. (2009). The 
CoRoT lightcurve covers almost 55 days from JD 2454138.1 to JD 2454192.8 
with 512 sec time sampling. We have removed data points with large error-bars 
which were typically upward outliers caused by high-energy particles hitting the 
CCD detector inside the measurement aperture. Most of these events occur 
when the satellite crosses the South Atlantic Anomaly (SAA). Elimination of 
these data points introduces small gaps in the lightcurve, with some of these 
gaps occurring quasi-periodically (due to SAA crossing). However, the number 
of discarded data points is relatively small (1168 out of 9229, 12.6%). To char- 
acterize the effect of periodic data rejection on the frequency analysis we have 
constructed the spectral window plot presented on Fig. 1. The amplitudes of 
the alias frequencies are of the order of 10 % or below. We have also corrected 
the lightcurve for a small (0.211 mmag/day) downward trend which is probably 
a result of gradual drift of the satellite pointing or some kind of instrumental 
decay. A detailed discussion of noise properties of the CoRoT data may be 
found in Aigrain et al. (2009) and Auvergne et al. (2009). 

Ground-based multicolor data 

In the framework of ground support to the CoRoT mission, we obtained Stromgren 
uvhy photometry of this object using the Wide Field Camera (WFC) on the 
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Figure 1: Spectral windows for the outlier-rejected data (light) and for the 
data excluding the primary eclipses (dark), which were used for frequency 
analysis. The inset shows a zoom-in to the lower frequencies. 



2.5 m Isaac Newton Telescope (INT) Roque Muchachos Observatory, La Palma, 
Canary Islands. The observations were conducted on JD 2454170.391 which 
corresponds to the orbital phase 0.197 of the binary system. The object was 
also observed in Johnson B and V filters by Deleuil et al. (2006) using the same 
telescope and camera. Infrared JHKg photometry of this object is available 
from the 2MASS catalog (Skrutskie et al. 2006). The 2MASS observations 
were conducted on JD 2451526.704 which corresponds to the orbital phase 
0.191 (by a lucky coincidence it is very close to the phase of our Stromgren 
photometry). The great advantage of the 2MASS photometry is that it was 
obtained using a camera which utilizes a beam splitter producing truly simulta- 
neous measurements in three filters. Therefore, 2MASS colors are not distorted 
by stellar variability. 

All available color information is summarized in Table 1. 



Interstellar extinction estimation 

In the absence of spectroscopic data it is important to constrain the influence 
of interstellar reddening, because information about intrinsic colors is necessary 
for unambiguous classification of the pulsating star. 

According to the standard tables by Schlegel et al. (1998) the total Galactic 
extinction in the direction of the object is large: Ay = 3™860, E{B — V) — 
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Table 1: Multicolor photometry of CoRoT 102980178 



Parameter 


Value 


Error 


Origin 


b — y = 


0.598 


±0.032 


INT/WFC 


y = 


15.723 


±0.034 


INT/WFC 


mi = 


-0.042 


±0.11 


INT/WFC 


Cl = 


0.891 


±0.086 


INT/WFC 


B = 


16.58 


±0.45 


INT/WFC 


V = 


15.73 


±0.19 


INT/WFC 


J = 


13.411 


±0.028 


2MASS 


H = 


12.759 


±0.023 


2MASS 


Ks = 


12.528 


±0.027 


2MASS 



1™164. The Schlegel et al. estimation is based on direct observations of far- 
infrared emission of interstellar dust, the material which causes absorption in 
visible and near-infrared bands. However, as it is noted by the authors, the 
extinction estimation for this position may be unreliable because of the low 
Galactic latitude (6 = — 1?50). There is an alternative way to estimate the 
Galactic extinction which is based on the relation between HI column density 
[Nhi) and extinction caused by dust: Nhi/Av = 1.79 x 10^^ cm^^ mag~-^ 
(Predehl &i Schmitt 1995). We estimate Nhi in the direction of the object using 
21 cm radio observations from the Leiden/Argentine/Bonn Galactic HI Survey 
(Kalberia et al. 2005, see also http://www.astro.uni-bonn.de/~webaiub/ 
english/tools_labsurvey .php). From these data we obtain A^^f/ = 0.682x 
10^^ cm^^ corresponding to Av{HI) — 3™810 which is in good agreement with 
the estimation obtained using the Schlegel et al. tables. 

Unfortunately, the distance to the binary system is unknown, as well as the 
exact distribution of the absorbing material along the line of sight. The intrinsic 
color of the system lies somewhere between the observed color and the color 
corrected for the total Galactic reddening along the line of sight. However, even 
reliable upper limits on the Galactic reddening will be useful for the following 
discussion. 

Interpretation and Modeling 

The lightcurve of CoRoT 102980178 shows eleven 1™26 deep primary and the 
equal number of 0^25 deep secondary eclipses (at phase 0.50). The nearly 
flat primary minimum bottom indicates that the eclipses are total. Ellipsoidal 
variability is evident in the out of eclipse lightcurve (Fig. 2). The Maximum 
preceding the primary minimum is 0^04 fainter then the maximum which follows 
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Figure 2: Detrended lightcurve of CoRoT 102980178 phased with hght ele- 
ments HJDuini = 2454139.0680+5^0548 x E. The solid line represents the 
F-primary model lightcurve. 

it (positive O'Connell effect, see Davidge &i Milone 1984 and Liu &i Yang 2003 
for a discussion of the effect). The period analysis with the Lafler & Kinnnan 
(1965) method leads to the following light elements: 

HJDmm = 2454139.0680 + 5^0548 x E 
±0.0007 ±0^0170 

Individual minima times were estimated using Kwee & van Woerden (1956) 
method and combined to produce a single primary minimum epoch. Its la- 
uncertainty was estimated from scatter of individual measurements on the Ob- 
served minus Calculated {O — C) plot produced with the above light elements. 
Period uncertainty was estimated following Schwarzenberg-Czerny (1991). 

Superimposed on this classical Algol-type lightcurve are OV^Ol oscillations 
with the period of 0^36372 (see Fig. 3 and the detailed discussion below). 
The oscillations are evident at all phases except during the primary minimum 
(Fig. 4). 

Ruling out an optical blend 

When analyzing a variable star showing two different types of variability (e.g. 
eclipses and pulsations) it is possible that the observed variability comes from 
two unrelated objects which just happened to be on the same line of sight. For 
the variable star described here, this possibility may be ruled out due to the 
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Figure 3: The residual lightcurve (F-primary model). 




Figure 4: The residual lightcurve (F-primary model) folded with the binary 
orbital period. The lower curve represents the section of residual lightcurve 
used for frequency analysis, it was shifted 0™1 along the vertical axis for 
visibility. 
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(a) CoRoT image 



(b) DSS finding chart 



Figure 5: (a) CoRoT image of the variable star, marked area is the measure- 
ment aperture. The image is in detector coordinates, east is up, north is to 
the left, (b) 5' x 5' DSS2 red image centered on CoRoT 102980178 (north 
is up, east is to the left). 



following reasons: 1) inspection of both the original CoRoT image (Fig. 5a) 
and the DSS image (Fig. 5b) reveal no other stars in the immediate vicinity 
of the variable and 2) pulsations disappear during the primary eclipse, which 
proves that the pulsating star is actually blocked from view. 

Oscillating component classification 

Following Moon (1986) we calculate reddening-free indices 

[ci] = ci - 0.19(5 - y) = 0.777 ± 0.087 

[toi] = TOi + 0.33(6 -y) = 0.155 ± 0.111 

which are consistent with a main sequence star of spectral class B8-9 or A7-F1 
within the uncertainty (see Fig. 6). 

The Johnson B and V band observations with INT/WFC by Deleuil et al. 
(2006) (see Table. 1) have large error bars and where not used in the analysis, 
however, we note that the {B — V) color of the star do not contradict our 
proposed spectral classification after accounting for the uncertainties in the 
Galactic extinction. 
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2MASS infrared colors can also be used to estimate spectral class. Account- 
ing for the uncertainty in the Galactic reddening discussed above, the infrared 
colors correspond to spectral class in the range early A - early M (Bessell &l 
Brett 1988). However, it is possible that the secondary component of the bi- 
nary system dominates the infrared light, so this spectral type estimation may 
correspond to the secondary star instead of the pulsating component. 

The multicolor photometry excludes early B spectral class, but is consistent 
with a reddened late B or F class star. 

The period and amplitude of oscillations are typical for 7 Doradus type 
but also within the possible range for young Slowly Pulsating B (SPB) and 
(3 Cephei stars. Information about the spectral class is crucial to distinguish 
between these possibilities: /3 Cephei variables are early-type B stars (BO - 
B2.5), SPB stars have spectral types between B2 and B9 (Stankov &i Handler 
2005) while 7 Doradus variables are typically early F type stars (e.g. Rodrfguez 
2002, Henry et al. 2007). f3 Cephei classification can be excluded on the basis 
of multicolor photometry. Note also, that according to Pigulski &l Pojmahski 
2009, f3 Cephei variables are characterized by periods < 0'?35 while the GCVS 
(Samus et al. 2009) gives 0^6 as the maximum period for this variability class. 
SPB classification is consistent with the Stromgren photometry, but SPB stars 
are usually characterized by longer pulsation periods of 0'?5 - 5*^ (Thoul 2009) 
and therefore, we consider classification of the oscillating component as an SPB 
star to be less likely. 7 Doradus stars are characterized by periods ranging from 
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0^4 to S'^ (Kaye et al 1999). Therefore, 7 Doradus type variable is the most 
probable classification of the oscillating component. 

Binary system modeling 

For the construction of the binary system model we have cleaned the lightcurve 
from outlier points and corrected it for the linear systematic trend. No prewhiten- 
ing for the pulsations was done. The oscillations should not affect the binary 
model because they are significantly smaller that the main features of the eclips- 
ing lightcurve. 

The light curve was analyzed with the PHOEBE software (Prsa Si Zwit- 
ter 2005, see also http://phoebe.fiz.uni-lj .si ), which provides a user- 
friendly interface and additional capabilities to the Wilson &i Devinney (1971) 
code. The software was recently modified to take into account CoRoT trans- 
mission functions, both for flux and limb darkening computations (for details 
see Maceroni et al. 2009). Being the orbit-unrelated variability of much smaller 
amplitude and shorter period than the eclipses, we just fitted the light curve 
folded in phase according to the ephemeris presented above (Fig. 2). 

After an initial screening of the parameter space, we proceeded to the light 
curve solution, using - at different stages - both the differential correction and 
the Nelder and Mead's Simplex algorithm. The fit adjustable parameters were: 
the inclination, i, the mass ratio, q = M2/M1, the surface potentials ^1^2 
(which, together with q, determine the components fractional radii ri 2), and 
the secondary effective temperature, Tes,2- The albedo Oj and gravity dark- 
ening coefficients /3i were fixed at the theoretical values (oi_2 = 1.0,0.5 and 
= 1.0,0.32, respectively, for stars with radiative or convective envelopes). 
We adopted a non-linear limb darkening law, whose coefficients were interpo- 
lated from the PHOEBE tables on-the-fly, according to the star temperature 
and surface gravity (logg). For the F-primary model (see below) we chose a 
logarithmic limb darkening law, for the B-primary one we chose a square-root 
law. The primary luminosity in the CoRoT passband, Li, was computed at 
each iteration in PHOEBE, rather than adjusted, to enhance convergence {L2 
is derived from the other parameters and model atmospheres). 

We started from a detached configuration but in all cases the iterative 
solution evolved towards a semi-detached configuration (with the less massive 
star in contact). We checked, as well, the effect of non-zero eccentricity, which 
however did not provide any improvement of the fit. 

It is well known that the solution of a single passband light curve is mainly 
sensitive to the Tes relative values, therefore a shift in Tes,i yields a shift in 
the same direction of Tes^2- It was therefore possible to fit the light curve in- 
corporating different assumptions about the temperature of primary (pulsating) 
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component: T^ff i = 7000 K (early F spectral class) and T^ff i = 11400 K (late 
B class). As expected, we obtained different temperatures of the secondary: 
Teff 2 = 4692 (K class) and Teff 2 = 6162 (late F) for the first and second model 
respectively. Apart from that, the models are very similar. The uncertainty of 
Teff 2 is dominated by our poor knowledge of Tgff 1. For a particular assumed 
value of Teff 1 the formal error bars of Teff 2 calculated from the least square fit 
are less than 10 K. 

To achieve an acceptable light curve solution the non-negligible O'Connell 
effect shall also be modeled, and for this purpose we introduced a bright spot on 
the primary component. Its parameters (location, size and temperature contrast 
factor, Tspot i/Teff 1) where not adjusted, but determined by trial and error, by 
putting the spot on the stellar equator and varying only the spot longitude and 
contrast factor (as the size and the latitude are highly correlated with them, see 
for example the discussion by Maceroni &; van't Veer 1993). It is reasonable 
to expect the existence of such a spot as the result of mass accretion from the 
secondary component. However different solutions (such as a dark spot on the 
colder secondary star, due in this case to surface activity), cannot be ruled out. 
The model parameters are summarized in Table 2. 

Table 2: Model parameters of the binary system 



Parameter Model with F primary Model with B primary 



assumed Teg 1 


7000 K 


11400 K 


Tefl 2 


4692 K 


6162 K 


q = M2/M1 


0.206 ±0.001 


0.184 ±0.001 


i 


84?61 ± 0.05 


85?61 ±0.04 


n 


0.1068 ± 0.0002 


0.1064 ± 0.0002 


r-2 


0.2447 ± 0.0003 


0.2363 ± 0.0003 


L2/L1 


0.7070 ± 0.0005 


0.7183 ±0.005 


^spot l/Teff 1 


1.22 


1.50 



The spot angular radius was fixed to 15° its colatitude to 90, the longitude 
of 240° and the temperature contrast factors in the table were found by 
trial and error. The uncertainties are formal fit errors (Icr). 

One of the lightcurve features which are not reproduced by the models is a 
small difference in depth among primary eclipses. The eleven observed eclipses 
are not sufficient to draw a conclusion about possible periodicity of the primary 
eclipse depth variation. The models also predict symmetric lightcurve shape 
during the total eclipse phase which is not observed. These two effects may be 
caused by starspots on the secondary component. 

As can be seen from the Fig. 2, there are some residual systematic deviations 
of the model with respect to observations. This will imply the presence of the 
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orbital frequency (forb) and its multiples (typically the even ones) in the residual 
lightcurve spectrum which is discussed in the following section. 

Pulsation frequency analysis 

The residual lightcurve resulting from subtraction of the binary model (version 
with early-F type primary) is presented on Fig. 3. It shows clear pulsation with 
a dominant frequency ^ 2.75 c/d. If folded with the binary orbital period (see 
Fig. 4) the residual lightcurve shows a rather complicated pattern instead of 
uncorrelated data. This means that the main pulsation frequency is close to an 
integer multiple of the binary orbital frequency (fl/forb = 13.898 ^ 14). The 
pulsation disappears during the primary eclipse (phases —0.05 - 0.05 on Fig. 4), 
and this phase interval was excluded from the frequency analysis. In order to 
check for possible artifacts introduced due to omitting parts of the lightcurve, 
the analysis was repeated including observations during total eclipses. The 
introduced artifacts were found to be negligible. 

The frequency analysis was performed using the SigSpec software (Reegen 
2007), leading to 99 formally significant frequencies including the CoRoT orbit 
and its harmonics. The first few frequencies are listed in the Table 3. The 
search for combination frequencies was performed using the software COMBINE 
developed by P. Reegen (see http://www.SigSpec.org). The frequency errors 
were estimated following Kallinger et al. 2008. The orbital frequency of the 
binary system was added to the frequency list, because from the plot of the 
resulting frequencies (Fig. 7) it is obvious that there are many combinations 
with the orbital frequency. 

Only two frequencies (fl and f3 in Table 3) appear to be genuine, which 
means not explained by a combination; f3 is close to forb but the frequency errors 
are small enough to conclude that the frequencies are significantly different (at 
lOfT level). 

Analysis of the residual lightcurve after subtraction of the B primary model 
leads to similar results (Table 3). The one major difference is that f8 becomes 
second most powerful harmonic in the spectrum. Its amplitude is nine times 
larger compared to residuals from the binary model with early-F type primary. 
fS in the B primary model residual spectrum can be associated with the orbital 
frequency (forb = 0.1978±0.0007), the difference between them is less then 4(t. 
Most probably, it is not a real frequency but an artifact introduced by imperfect 
binary model subtraction. 

All the other frequencies can be explained by combinations mainly with the 
orbital frequency. This is what can be expected from a close binary because 
the primary component is likely to be distorted by tidal forces. If the F type 
primary hypothesis is accepted, the solution for the independent frequencies has 
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Table 3: Frequencies detected in the residual lightcurve 





Model with F 


primary 


Ti JT 11 ■ J 1 T> 

Model with B 


primary 




Frequency 


Amplitude 


Frequency 


Amplitude 




[c/d] ' 


[mmag] 


[c/d] 


[mmag] 


fl 


2.7494 ±0.0008 


9.52 ±0.42 


2.7494 ± 0.0009 


9.68 ± 0.48 


(4forb) f2 


0.7917 ±0.0008 


7.48 ± 0.35 


0.7920 ± 0.0009 


6.31 ±0.32 


f3 


0.2115 ±0.0012 


3.28 ±0.25 


0.2181 ±0.0021 


1.97 ±0.23 


(2forb) f4 


0.4004 ± 0.0015 


3.19 ±0.26 


0.4012 ± 0.0016 


2.57 ±0.22 


(lOforb) f5 


1.9784 ±0.0017 


3.27 ±0.30 


1.9783 ±0.0014 


3.86 ±0.30 


(8forb) f6 


1.5823 ±0.0017 


2.81 ±0.27 


1.5822 ±0.0013 


4.17 ±0.29 


(isforb) fr 


2.5657 ±0.0019 


2.48 ±0.25 


2.5650 ±0.0019 


2.32 ±0.25 


(forb) f8 


0.1978 ± 0.0031 


1.00 ±0.17 


0.1940 ±0.0008 


9.05 ± 0.40 



two components: fl and fS. 

Conclusions 

We have identified CoRoT 102980178 as a semidetached eclipsing binary with 
a 5^0548 orbital period. It shows an Algol-type lightcurve with 1™26 deep 
primary and 0™25 deep secondary eclipses (in CoRoT photometric band). 

Our lightcurve modeling leads to a semidetached configuration with the 
photometric mass ratio q = 0.2 and orbital inclination i = 85°. We note, that 
the inverse problem of binary model fitting using only a broad passband light 
curve is degenerate and the proposed solution might not be unique. 

The more massive component of the system is pulsating with the primary 
frequency of 2.7494 c/d (0':'36372 period) and 0™01 amplitude. The observed 
period and amplitude of pulsations as well as ground-based multicolor photom- 
etry favor its classification as 7 Doradus type variable, however, based on our 
data we cannot exclude the SPB-hypothesis for the primary. The detailed fre- 
quency analysis suggests the presence of an additional independent pulsation 
mode with the frequency of 0.21-0.22 c/d. Other detected frequencies can be 
explained as combinations of the above frequencies and the orbital frequency. 

The system deserves a detailed spectroscopic study which could better con- 
strain physical parameters of the pulsating component. A Stromgren photom- 
etry of the primary eclipse could also be useful for this purpose. 
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